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ABSTRACT 

We have investigated the gas-phase chemistry in dense cores where ice mantles containing ethanol and 
other alcohols have been evaporated. Model calculations show that methanol, ethanol, propanol, and butanol 
drive a chemistry leading to the formation of several large ethers and esters. Of these molecules, methyl ethyl 
ether (CH 3 OC 2 H 5 ) and diethyl ether [(C 2 H 5 ) 2 0] attain the highest abundances and should be present in 
detectable quantities within cores rich in ethanol and methanol. Gas-phase reactions act to destroy evaporated 
ethanol and a low observed abundance of gas-phase C 2 H 5 OH does not rule out a high solid-phase abun- 
dance. Grain surface formation mechanisms and other possible gas-phase reactions driven by alcohols are 
discussed, as are observing strategies for the detection of these large interstellar molecules. 

Subject headings: ISM: abundances — ISM: clouds — ISM: molecules — molecular processes 


1. INTRODUCTION 

Methanol (CH 3 OH), the simplest alcohol, is widely observed 
in both the gas phase and the solid phase of the dense inter- 
stellar medium (e.g., Irvine, Goldsmith, & Hjalmarson 1987; 
Tielens & Allamandola 1987). Infrared observations show that 
it can be a significant fraction of molecular ice mantles relative 
to water (Tielens & Allamandola 1987). Gas-phase abundances 
relative to H 2 range from ~ 10 9 in cold, dark clouds (Friberg 
et al. 1988) to 10” 7 to 10” 6 in the hot molecular cores associ- 
ated with regions of massive star formation (Menten et al. 
1988) and ultracompact H II regions (Churchwell 1990). It is 
now accepted that the highest gas-phase methanol abundances 
are found in regions where ice mantles have recently evapo- 
rated from warm dust grains. Gas-phase processes cannot 
produce the observe abundances of 10“ 7 or more (Millar, 
Herbst, & Charnley 1991a) and solid state chemical reactions 
on grains are necessary. 

Ethanol (C 2 H 5 OH), second in the homologous series, is not 
observed in dark clouds (Irvine et al. 1987) but it is present in 
star-forming regions with gas-phase abundances estimated to 
lie in the range of ~ 10” 8 (Millar et al. 1988) to 10 6 (Turner 
1991). The lowest abundances cannot be reproduced in gas- 
phase chemistry models (Herbst & Leung 1989; Millar et al. 
1991a), even under the most favorable conditions (Charnley, 
Tielens, & Millar 1992). Observational spectra of the C-H 
stretching band near 3.4 fxm in the diffuse interstellar medium 
bear a striking resemblance to that of solid n-butanol (n-butyl 
alcohol) at 10 K (Sandford et al. 1991). Thus, the observational 
evidence is that interstellar CH 3 OH and C 2 H 5 OH can be 
formed by grain surface chemistry much more efficiently than 
by gas-phase chemistry, and that even higher alcohols could be 
present in grain mantles. 

Recent theoretical work has demonstrated that the chem- 
istry of hot molecular cores is active (Charnley & Tielens 1992; 
Charnley et al. 1992; Caselli, Hasegawa, & Herbst 1993; 
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Charnley & Millar 1994): evaporated ices drive selected gas- 
phase pathways depending upon the relative molecular com- 
position of the ices. This scenario can explain the hot core 
abundances of complex molecules such as dimethyl ether 
QCH 3 ) 2 0], methyl formate (HCOOCH 3 ), ketene (CH 2 CO) 
and acetaldehyde (CH 3 CHO), as well as the molecular differ- 
entiation between individual cores. Dimethyl ether and methyl 
formate are not the products of grain surface chemistry but are 
formed by ion-molecule reactions from evaporated methanol; 
the former by the well-known self-alkylation process (e.g., 
Karpas & Mautner 1989) and the latter by alkyl transfer to 
formaldehyde from protonated methanol. These observational, 
theoretical, and experimental considerations lead one to 
inquire as to the nature of astrochemistry in sources where 
large quantities of alcohols are present. 

In this paper we adopt the view of Charnley et al. (1992) that 
only specific classes of molecule, in particular the alcohols, are 
produced in grain mantles and evaporated in hot core regions. 
In this way, we are able to produce quantitative models of the 
subsequent chemistry and suggest particular observational 
tests of the theory. In § 2 we outline the theoretical model and 
describe an ion-molecule chemistry, based on experimental 
measurements involving methanol, ethanol, propanol, butanol, 
and formaldehyde, as well as the complex molecules that may 
be formed from them. In § 3, we present model results of calcu- 
lations of hot core chemistry under conditions appropriate to 
the Compact Ridge component of Orion A. In § 4, we discuss 
the viability of our assumptions concerning alcohol formation 
on grains, other possible routes to gas-phase complexity 
attainable via alcohol chemistry, and briefly outline prospects 
for observing large organic molecules. Our conclusions are 
summarized in § 5. 

2. THEORY 

2.1. Model 

We have employed a chemical kinetic model to study hot 
core chemistry (see Brown, Charnley, & Millar 1988) when ices 
are evaporated containing significant quantities of ethanol, iso- 
propyl alcohol (/-C 3 H 7 OH), «-propyl alcohol («-C 3 H 7 OH), 
t-butyl alcohol (f-C 4 H 9 OH), and n-butyl alcohol (n-C 4 H 9 OH), 
as well as methanol and formaldehyde. We adopted the 
approach taken by Charnley et al. (1992) to calculate the 
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chemical evolution in the hot static gas. We consider an initial 
state where all the heavy atomic and molecular material is 
incorporated in mantle molecules and where the gas phase 
consists of molecular hydrogen and helium atoms along with 
their related molecular and atomic ions. The molecular ice 
mantles are composed of H 2 0, CO, N 2 , 0 2 , CH 4 , C 2 H 4 , 
C 2 H 6 , H 2 CO, CH 3 OH, C 2 H 5 OH, i-C 3 H 7 OH, n-C 3 H 7 OH, 
f-C 4 H 9 OH, and n-C 4 H,OH. Evaporation of these ice mantles 
drives a hot gas-phase chemistry and the evolution of molecu- 
lar complexity. 

We consider a dense core with n H = 2 x 10 cm and tem- 
perature of 100 K, conditions appropriate for the Orion 
Compact Ridge (Blake et al. 1987). The mantle abundances 
suggested by Charnley et al. (1992) yield good agreement for 
the Compact Ridge abundances of several organic molecules 
and we take this source as prototypical of other methanol-rich 
hot cores. For the calculations of this paper, the input mantle 
abundances for the first nine molecules listed above were from 
Table 1 of Charnley et al. (1992); that is, 2 x 10 s ,lx 10 , 
7 x 1(T\ 2 x 10~ 7 , 1 x 10 7 , 8 x 10 8 , 1 x 10 , 

1.2 x 10 7 , and 1 x 10 ” 6 , respectively. These values were fixed 
in all models. To span the range of observational estimates, we 
performed three calculations with input ethanol abundances of 
1 x 10 6 , 1 x 10 7 , and 1 x 10' 8 , and present results for the 
first two. Note that having no ethanol in the ices leads to 
ethanol gas-phase abundances of less than 10 11 (Charnley et 
al. 1992). In each of the three models the total input propanol 
and butanol abundances were taken to be respectively 10 and 
100 times less than that of ethanol, divided equally between the 
two isomeric forms considered here. The adopted values are 
based upon the reasonable (but unverifiable) assumption that 
the mantle abundances of higher homologues fall off with 
increasing size. 

2.2. Interstellar Alcohol Chemistry 

The basic gas-phase chemistry employed is well documented 
(Millar et al. 1991a, b; Charnley et al. 1992). For the chemistry 
of the higher alcohols we have used the experimental results of 
Karpas & Mautner (1989). The basic reaction process is alkyl 
transfer from a protonated alcohol to a neutral one leading to 
the protonated pure or mixed ether (i.e., containing identical or 


different alkyl groups). The principal ion-molecule reactions 
involving interstellar alcohols and the rate coefficients adopted 
are listed in Table 1. Figure 1 shows schematically the main 
reactions of interest in the coupled methanol-ethanol chem- 
istries. Each of HCO + , H 3 0 + , and CH 3 OH 2 + can protonate a 
higher neutral alcohol but only CH 3 OH is protonated by H 3 . 
Ethanol is destroyed by H 3 (Millar et al. 1991b) and we 
assume, pessimistically, that this is also the case for higher 
alcohols; it is certainly true for dimethyl ether and diethyl ether 
(Lee, Drucker, & Adams 1992). The reactions leading from 
n-C 3 H 7 OH and n-C 4 H 9 OH to the pure ethers, (n-C 3 H 7 ) 2 0 and 
(n-C 4 H 9 ) 2 0, are analogous to those depicted for CH 3 OH in 
Figure 1, except that their protonated forms do not proton 
transfer to other alcohols. The normal isomers of propanol and 
butanol apparently do not form mixed ethers (Karpas & 
Mautner 1989); however, in the case of the isomers i-C 3 H 7 OH 
and r-C 4 H 9 OH, their protonated forms can react with 
CH 3 OH leading to the production of mixed ethers. The chem- 
istries of i-C 3 H 7 OH, and t-C 4 H 9 OH are each analogous to 
that depicted for ethanol in Figure 1. 

Karpas & Mautner (1989) showed that as the degree of 
alkylation (i.e., the total number of alkyl groups in the reacting 
system) increased, the rate coefficients of alkyl transfer reac- 
tions tended to decrease. We have thus ignored many reactions 
which contain three or more alkyl groups such as those 
between protonated ethers and alcohols. Of these, two reac- 
tions that could be important are 

(CH 3 ) 2 OH + + CH 3 OH -► (CH 3 ) 3 0 + + H 2 0 , (1) 

(CH 3 ) 2 OH + + C 2 H 5 OH CH 3 OC 2 H 6 + + CH 3 OH (2) 

Reaction (2) is endothermic (Karpas & Mautner 1989) and 
proton transfer is preferred. Neglect of reactions like reaction 
(1) greatly reduced the number of complex species added to the 
original reaction network. We have, however, included three 
reactions involving three alkyl groups: reactions of protonated 
mixed ethers with the neutral alcohols from which they were 
derived (reactions 4, 5, & 6 of Table 1). These reactions lead to 
the pure ether and are potentially important in determining the 
relative abundances of pure and mixed ethers. Note that 
CH 3 OH 2 does not react with C 2 H 5 OH or higher alcohols to 



Fig. 1 —Schematic of the principal ion-molecule reaction paths leading from methanol and ethanol to ethers 
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TABLE 1 


Ion-Molecule Reactions Involving Methanol, Ethanol, Propanol, and Butanol 


Reaction 


1 

C 2 H 5 OH 2 + 

+ h 2 co 

-* HCOOC 2 H 6 + 

+ h 2 

2 

C 2 H,OH 2 + 

+ ch 3 oh 

- CH 3 OC 2 H 6 ‘ 

+ h 2 o 

3 

c 2 h 5 oh 2 + 

+ C 2 H 5 OH 

-(C 2 H,) 2 OH* 

+ h 2 o 

4 

CHjOC 2 H* 

+ C 2 H 5 OH 

->(C 2 H,) 2 OH + 

+ CH30H 

5 

i-c 3 h 7 och; 

+ i-C 3 H 7 OH 

-(i-C 3 H 7 ) 2 OH + 

+ CH30H 

6 

r-C 4 H 9 OCH; 

+ r-C 4 H 9 OH 

-»(r-C 4 H,) 2 OH + 

+ CH30H 

7 

h-c 3 h,oh; 

+ h 2 co 

-► »-C 3 H 7 OOCH 2 + 

+ H 2 

8 

n-C 3 H 7 OH 2 + 

+ m-C 3 H 7 OH 

(n-C 3 H 7 ) 2 OH + 

+ h 2 o 

9 

j-C 3 H 7 OH 2 + 

+ h 2 co 

- i-C 3 H 7 OOCH 2 + 

+ h 2 

10 

i-CjHjOH,* 

+ CH30H 

-m-c 3 h 7 och; 

+ h 2 o 

11 

i-C 3 H 7 OH 2 + 

+ i-C 3 H 7 OH 

-0C 3 h 7 ) 2 oh + 

+ h 2 o 

12 

n-C t H,OH; 

+ rt-C 4 H 9 OH 

-+ (n-C 4 H 9 ) 2 OH + 

+ h 2 o 

13 

n-C,H,OHj* 

+ h 2 co 

-► n-C 4 H 9 OOCH 2 + 

+ h 2 

14 

t-C 4 H,OH 2 

+ ch 3 oh 

- £-C 4 H 9 OCH 4 

+ h 2 o 

15 

f-C 4 H,OH 2 

+ h 2 co 

- f-C 4 H 9 OOCH l 

+ h 2 

16 

(-C.H,OH 2 + 

+ £-C 4 H 9 OH 

-* (f-C 4 H 9 ) 2 OH + 

+ h 2 o 

17 

CHjOH 2 + 

+ c 2 h 5 oh 

->C 2 H,OH 2 * 

+ CH30H 

18 

ch 3 oh 2 + 

+ c 2 h,oh 

-►c 2 h 5 och; 

+ h 2 o 

19 

ch 3 oh 2 * 

+ i-C 3 H 7 OH 

-► i-C 3 H 7 OH 2 

+CH30H 

20 

CHjOH 2 + 

+ n-C 3 H 7 OH 

n-C 3 H 7 OH 2 

+ ch 3 oh 

21 

ch 3 oh 2 + 

+ n-C 4 H 9 OH 

- «-C 4 H 9 OH 2 

+ CH30H 

22 

ch 3 oh 2 + 

4- r-C 4 H 9 OH 

- f-C 4 H 9 OH 2 

+ CH30H 


NoTE:a(~b) -ax 10 “*. 

References. — ( 1 ) this work ; (2) Karpas & Mautner 1 989. 

b R eac t| on observed by Karpas & Mautner; their estimate of k, our product distribution. 

Reaction observed by Karpas & Mautner; our estimate of k. 
c Reaction postulated by Karpas & Mautner; our estimate of k. 


k(cm 3 s *) 


1 .70( — 11) 
3.00( — 1 1) 
6.80( — 1 1) 
I.00(— 11) 
1.50(— 11) 
2.25( — 1 1) 
1.00(-1I) 
1 . 00 (- 10 ) 
4.30(-n) 
3.00( — 11) 
6 . 10 (— 10 ) 
1.50( — 10) 
1 . 00 (- 11 ) 
3.00( — 11) 
1 . 00 (- 11 ) 
6.60(— 10) 
3.00(— 11) 
0 . 00 ( + 00 ) 
3.00( — 11) 
3.00( — 11) 
3.00(— 1 1) 
3.00(— 1 1) 


Reference 


1 “ 

l b 

2 

V 

T 

2 b 

2 * 

2 

1“ 

l c 

2 

2 

1 

l c 

1 

2 

2 b 

2 

1 

1 

1 

1 


form protonated mixed ethers, and that instead proton transfer 
occurs. As CH 3 OH is the most abundant alcohol in hot cores, 
we have included only reactions between it and the protonated 
isomers of propanol and butanol. We have not considered any 
reactions between complex species derived from ethanol, 1 - 
C 3 H 7 OH, n-C 3 H 7 OH, t-C 4 H 9 OH, or n-C 4 H 9 OH 
Table 1 shows that rate coefficients for alkyl transfer 
increase with increasing alcohol size and degree of branching. 
Experiments (Karpas & Mautner 1989) and theory (Herbst & 
Dunbar 1991) show that the rate coefficients for radiative 
association increase with reactant size. We have ignored radi- 
ative association reactions since the actual abundances of the 
molecules in whose chemistry it may play a significant part, 
propanol and butanol, are unknown. Detailed consideration of 
such extremely complex molecules as might result from this 
route is beyond the scope of this study. Important, and tenta- 
tive, inclusions in Table 1 are that, in analogy with the reaction 
believed to lead to interstellar methyl formate (e g., Millar et al 
1991b), that is, 

CH 3 OH 2 + + H 2 CO - HCOOCH 4 + + H 2 , (3) 

we have assumed that reactions of protonated higher alcohols 
with formaldehyde, for example, 

C 2 H 5 OH 2 + + H 2 CO - HCOOC 2 H 6 + + H 2 , (4) 

lead to the protonated ester. Recombination of the molecular 
ion produced in reaction (4) should form ethyl formate 
(HCOOC 2 H 5 ). We have been unable to find experimental veri- 
fication of this assumption and note that Karpas & Mautner 
(1989) give that the products of reactions (3) and (4) include 
CH 3 OCH 2 and C 2 H 5 OCH 2 , respectively. Inclusion of these 
molecular ions and their chemistries increases the complexity 
of the reaction network and, as we are motivated to explore the 


possibility of producing large molecules by reactions analo- 
gous to reaction (3), we retain the reactions of Table 1 but note 
that our results may be compromised beforehand. 

Finally, to the reactions of Table 1 we have added appropri- 
ate destruction reactions for complex organic molecules. 
Protonation of large neutral molecules by H 3 , HCO + , and 
H 3 0 + proceeds with a rate coefficient of 2-3 x 10" 9 cm 3 s" *. 
Charge transfer by H + and C + occurs at a rate 3 x 10" 9 cm 3 
s' 1 . Destruction by He + , with a rate coefficient of 3 x 10‘ 9 
cm 3 s" ! , is assumed to cleave the large organic molecules to 
yield several simple hydrocarbons in the product distribution 
and so to limit the development of molecular complexity. Dis- 
sociative recombinations occur at a rate of 1.5 x 10" 7 cm 3 
s' 1 . For the protonated ethers and esters derived from ethanol, 
propanol, and butanol, equal branching ratios were assumed 
for the channel yielding the original alcohol and the channel 
producing the associated neutral ether or ester. The total 
number of reactions added to the original network was 171. 

3. RESULTS 

Understanding the chemistry driven by alcohols is straight- 
forward (Fig. 1). Cosmic-ray ionization leads to H 3 , HCO+, 
and H 3 CU in the core and these protonate the neutral alco- 
hols. The ion-molecule chemistry then proceeds by reactions 
between protonated and neutral alcohols (principally self- 
alkylation), leading to protonated pure and mixed ethers; dis- 
sociative recombination then produces the neutral ether 
molecule. As the core evolves the alcohols are gradually 
destroyed and the derived ethers and esters concomitantly dis- 
appear. 

Figure 2 shows the evolution of the chemistry driven by 
evaporated methanol and that of the principal molecular ions; 
the evolution is the same in each of the models described 
below. The calculated abundances of (CH 3 ) 2 Q and 
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Fig. 2. — Chemical evolution of the methanol-formaldehyde chemistry and 
the major molecular ions. The evolution is almost identical in all models as the 
input CH 3 OH and H 2 CO abundances are the same in each. The physical 
conditions were n H = 2 x 10 6 cm” 3 and T — 100 K. The mantles were injected 
into the gas at t = 0. See text for further details. 

HCOOCH3 can be used to monitor the evolution of complex- 
ity in hot cores; for example, in the Orion Compact Ridge 
these molecules appear to be peaking in abundance (Charnley 
et al. 1992), implying an age of 10 4 -10 5 yr. Figures 3 and 4 
show the results of two model calculations in which the mantle 
abundances of ethanol, propanol, and butanol differ by a 


factor of 10. As expected, the abundances of the derived species 
are proportional to the input abundance of the relevant 
alcohol. The chemistries of the higher alcohols evolve roughly 
on a similar timescale as the methanol chemistry although, 
since they are destroyed and not protonated by H3 , the higher 
alcohols are more rapidly removed than methanol. Apart from 
being the only isomeric forms that react with CH 3 OH 2 , i- 
C 3 H 7 OH and f-C 4 H 9 OH react faster with their protonated 
counterparts than do the corresponding n-isomers (Table 1). 
This results in higher abundances of species containing the 
isopropyl and tertiary butyl groups in the models. On the one 
hand, the abundance of methanol is much higher than that of 
the other alcohols in the models, enhancing the production 
rate of mixed over pure ethers. On the other hand, the reaction 
rate coefficients involved in the production of the pure ethers is 
larger than those for mixed ethers. These two opposing effects 
balance for the models shown in Figure 3 and so the abun- 
dances of diethyl ether [(C 2 H 5 ) 2 0] and methyl ethyl ether 
(CH 3 OC 2 H 5 ) are comparable. At lower input abundances for 
ethanol (Fig. 4) CH 3 OC 2 H 5 always dominates. A similar argu- 
ment holds for the pure and mixed ethers formed from the 
branched isomers of propanol and butanol. 

We now apply these theoretical results to the interpretation 
and molecular phenomenology of hot core sources. The 
maximum observed abundances of CH 3 OH and C 2 H 5 OH are 
determined by grain surface chemistry; given our basic 
premise, it is tautological to state that both propanol and 
butanol may also be present in hot cores. Considering the 
higher alcohols first, and taking a fractional abundance of 
10” 11 as a lower limit, then, unless the mantle abundances of 
propanol and butanol are significantly higher than we have 
assumed, only a few large ethers can conceivably be formed at 
detectable abundances. These are (r-C 3 H 7 ) 2 0 (di-isopropyl 
ether) and i-C 3 H 7 OCH 3 (methyl isopropyl ether); t- 
C 4 H 9 OCH 3 (methyl r-butyl ether) is the most abundant 



Fig. 3.— Chemical evolution of the (a) ethanol, (b) propanol, and (c) butanol chemistries in the hot core model depicted in Fig. 2. The input abundances of each 
alcohol were 1 x 10“ 6 , 1 x 10 -7 ,andl x 10 “ 8 , respectively. 
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Fig. 4. — Same as Fig. 3 exept that the input abundances of ethanol, propanol, and butanol were 1 x 10 7 , 1 x 10 _8 ,andl x 10 ” 9 , respectively 


complex molecule formed from butanol. We expect that these 
alcohols may play only a limited part in producing (detectable) 
molecular complexity in the gas phase. By contrast, for 
ethanol, a range of abundance estimates are available (~ 10“ 6 
to 10“ 8 ; Blake et al. 1987; Millar et al. 1988, 1995; Turner 
1991), and our calculations demonstrate that this molecule can 
produce appreciable abundances of (C 2 H 5 ) 2 0 and 
CH 3 OC 2 H 5 over most of this range. Of the esters considered 
in this study only HCOOC 2 H 5 (ethyl formate) attains a signifi- 
cant abundance, and this in the models containing the most 
ethanol (see Fig. 3). 

In the Compact Ridge the observed (CH 3 ) 2 0 abundance is 
about 10 8 (Blake et al. 1987), implying that the time since 
mantle evaporation lies in the range ~ 10 4 -10 5 yr (Fig. 2). If 
Turner’s high estimate of the ethanol abundance ( — 10 _ a ; 
Turner 1991) is correct, both CH 3 OC 2 H 5 and (C 2 H 5 ) 2 0 
should be present in this source with abundances of a few times 
10' 9 . If the ethanol abundance of Blake et al. (1987) is more 
appropriate (~10 -11 ), then the models predict abundances 
below detection limits. In W51 M, Millar et al. (1988) found 
that (CH 3 ) 2 0 has the same abundance as in the Compact 
Ridge but the C 2 H 5 OH abundance is 10 8 . Millar et al. inter- 
preted these observations with a model in which the products 
of grain surface chemistry are sensitive to the H-atom abun- 
dance. However, it is difficult to use this model to predict the 
ethanol/dimethyl ether ratio in particular sources, so we inter- 
pret the observations in terms of our detailed numerical calcu- 
lations. The transient nature of the ethanol-driven chemistry 
leads to two possible interpretations of the W51 M observ- 
ations. First, in a model which initially contains 100 times 
more ethanol than presently observed in W51 M (i.e., Fig. 3), 
(CH 3 ) 2 0 and C 2 H 5 OH attain observed abundances after 
about 10 5 yr and both CH 3 OC 2 H 5 and (C 2 H 5 ) 2 0 have abun- 
dances of ~2 x 10" 10 . The important point here is that 
ethanol may originally have been very abundant (in the evapo- 


rated mantles) but has been subsequently destroyed by gas- 
phase reactions. Despite much lower observed ethanol 
abundances, the daughter products may still be present. 
Clearly, sources with very large ethanol abundances should 
show the largest ether abundances. Second, if ethanol is 
injected from grains at its observed abundance in W51 M, then 
we calculate that the peak ether abundances are severely 
depressed. In particular, the peak CH 3 OC 2 H 5 abundance falls 
from ~4 x 10“ 9 (Fig. 3) to ^4 x 10 11 ; all species derived 
from the higher alcohols have abundances much lower than 
CH 3 OC 2 H 5 . Detection of either CH 3 OC 2 H 5 or (C 2 H 5 ) 2 0 in 
W51 M would confirm the importance of ethanol as a grain 
mantle constituent. 

4. DISCUSSION 

We now address the question of the formation of interstellar 
alcohols by solid state chemical reactions, the detectability of 
complex organic molecules, and other aspects of alcohol chem- 
istry. 

4.1. Formation of Alcohols on Grains 

We briefly discuss three possible reaction pathways, with 
particular reference to the synthesis of methanol and ethanol, 
where observation clearly demands such a process. 

First, we consider the importance of H and O atom attack 
on species containing carbon double and triple bonds. Such 
reactions generally have activation barriers but quantum tunel- 
ling or thermally promoted reactions are still possible on grain 
surfaces (Tielens & Hagen 1982; Tielens 1992). Hydrogen atom 
addition to CO may lead to CH 3 OH and matrix isolation 
studies show that HCO is readily formed at 10 K (van IJzen- 
doorn 1985; Bohn, Allamandola, & Sandford 1995). Forma- 
tion of higher alcohols by this route requires intermediate oxi- 
dation of alkyl radicals. For example, H addition toC 2 H 2 with 
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oxidation of either the vinyl or ethyl radicals could produce 
ethanol (Tielens 1992). Acetylene is predicted to be highly 
abundant (%10 -6 ) on grain surfaces within dense conden- 
sations by direct accretion from the gas phase (Brown et al. 
1988). Evans, Lacy, & Carr (1991) have observed large abun- 
dances of C 2 H 2 in the Orion hot core where its origin must be 
attributed to evaporation from grains (Charnley et al. 1992). 
Accretion of other acetylenic compounds, followed by hydro- 
genation and oxidation, are plausible surface routes to various 
higher alcohols. Hot core observations of methyl acetylene 
(CH 3 CCH) and diacetylene (C 4 H 2 ), possible precursors of 
propanol and butanol, are not available. Although CH 3 CCH 
is detected in a number of molecular clouds (Irvine et al. 1981 ; 
Askne et al. 1984; Andersson, Askne, & Hjalmarson 1984), and 
is sometimes seen toward hot core sources (e.g., G34.3), it is 
always found to have a low rotational temperature — there is 
no observational evidence that CH 3 CCH is a “hot core” mol- 
ecule (Ziurys & McGonagle 1993). Grain formation of ethanol 
and dimethyl ether was considered by Millar et al. (1988). It is 
now clear that (CH 3 ) 2 0 is formed by gas-phase chemistry. In 
the qualitative model devised by Millar et al. (1988), the pro- 
posed route to C 2 H 5 OH can be discounted since it relies on an 
unrealistically low abundance of H atoms in the accreting gas. 

Second, we consider radical-radical surface reactions. At 
dust temperatures of 10 K, only O, C, and N atoms can diffuse 
rapidly enough to participate in the surface chemistry. Rad- 
icals, such as OH and CH 3 , are expected to remain trapped in 
their residence site (Tielens & Hagen 1982) and hence are effi- 
ciently hydrogenated to water and methane, with methanol 
production very limited (see Tielens & Hagen 1982; Brown & 
Charnley 1990). At higher dust temperatures ( ~ 20 K), OH and 
CH 3 radicals can diffuse on a grain surface and their reaction 
may lead to methanol. This scenario predicts variations in the 
concentration of methanol in interstellar grain mantles 
depending on the accretion temperature of the grains (cold: 
water-rich; warm: methanol-rich). Accretion, oxidation, and 
hydrogenation of CH 3 would also be important in methanol 
production; formation of higher alcohols would require the 
oxidation of accreted ethyl, propyl, and butyl radicals and is 
therefore contingent on the existence of a well-developed 
hydrocarbon chemistry in the accreting gas. 

Third, we consider diffusion-controlled reactions in bulk 
molecular ices. The icy mantles of cold dust grains in molecular 
clouds contain the products of H atom addition surface reac- 
tions: H 2 0, NH 3 , CH 4 , H 2 S, as well as other molecules 
present by accretion or by hydrogenation of accreted mol- 
ecules: CO, N 2 , 0 2 , C 2 H 2 , C 2 H 4 , C 2 H 6 . Ultraviolet photoly- 
sis of these simple ices, for example by the Prasad-Tarafdar 
flux (Prasad & Tarafdar 1983), can lead to atoms and simple 
radicals being frozen in the amorphous ice matrix at dust tem- 
peratures of 10 K: OH, CH 3 , NH 2 , CN, C 2 H 3 , C 2 H 5 (see 
Hagen, Allamandola, & Greenberg 1979; Greenberg 1982). 
Warming of these ices in star-forming regions will promote 
radical diffusion, annealing, and the formation of more 
complex molecules (Agarwal et al. 1985; Schutte 1988). Reac- 
tions of the hydroxyl radical with methyl and ethyl radicals 
could produce methanol and ethanol. If, during formation of a 
cold dense core, H addition to C 2 , C 3 , and C 4 , along with 
other partially hydrogenated carbon chain species, is indeed 
efficient (Hasegawa, Herbst, & Leung 1992), then UV photoly- 
sis and warming of ices containing the resulting aliphatic 
hydrocarbons such as propane (C 3 H 8 ) and butane (C 4 H 10 ) 
could produce higher alcohols. 


4.2. Other Alcohol-driven Reactions 

This work has been concerned primarily with ion-molecule 
reactions involving simple alcohols and formaldehyde. It is 
natural to inquire what complex molecules could be driven by 
reactions of CH 3 OH 2 and C 2 H 5 OH 2 with other abundant 
mantle molecules such as NH 3 and HCN. Hydrogen cyanide 
can be formed by gas-phase reactions in hot cores (Charnley et 
al. 1992) but is also probably abundant on grains prior to 
evaporation (Brown et al. 1988), as witnessed by the presence 
of DCN in hot cores (Brown & Millar 1989; Mangum, Plam- 
beck, & Wootten 1991). Two reactions of interest are 

CH 3 OH 2 + + HCN - CH 3 NCH + + H 2 0 , (5) 

C 2 H 5 OH 2 + + HCN - C 2 H 5 NCH+ + H 2 0 . (6) 

Both reactions (5) and (6) have been studied in the laboratory 
and have measured rate coefficients of 1.8 x 10“ 11 cm 3 s _1 
and 1.5 x 10“ 11 cm 3 s’ 1 , respectively (Mautner & Karpas 
1986). Two interesting possibilities exist, depending upon 
whether reactions (5) and (6) lead to isocyanides or to cyanides. 
Mautner & Karpas (1986) have argued that protonated iso- 
cyanides are formed and hence we would expect, depending on 
the branching ratios for dissociative recombination, that 
CH 3 NC and C 2 H 5 NC would be viable species for detection in 
hot cores. In this context, we note that CH 3 NC has been tenta- 
tively detected in TMC-1 by Irvine & Schloerb (1984) and 
claimed in Sgr B2 by Cernicharo et al. (1988). An alternative 
possibility is that, following formation of the protonated iso- 
cyanides, isocyanide-cyanide isomerization occurs (see Defrees, 
McLean, & Herbst 1985) and that the ultimate neutral pro- 
ducts are CH 3 CN and C 2 H 5 CN. More CH 3 NC measurements 
of methanol-rich cores are needed. The nondetection of iso- 
cyanides in hot cores would support the second possibility. 
Alkyl cation transfer reactions may in fact account for most of 
the organic molecules seen in hot cores (Charnley 1995). For 
systems involving a large number of atoms, radiative associ- 
ation can become significant in the formation of large molecu- 
lar ions (Herbst & Dunbar 1991), leading to increased 
molecular complexity. However, the rate coefficients involved 
in radiative association are usually a strong function of tem- 
perature, indicating that, to a first approximation, such pro- 
cesses can be neglected in hot gas chemistry. 

4.3. Detection Prospects for Large Organic Molecules 

It is of some interest to consider the possibility of detecting 
the large molecules predicted here to be present in hot cores. 
Microwave studies of large molecules are difficult. For large 
species such as these, centrifugal distortion is important and 
hence transition frequencies cannot be easily predicted from 
rotational constants but have to be measured directly. The 
small rotational spacing leads to large partition functions, par- 
ticularly at hot core temperatures, and the intensity in any line 
is small, more so for nonlinear molecules. However, we can get 
an indication of the expected antenna temperatures of some 
large molecules by noting that, for optically thin emission, the 
intensity of a given transition is proportional to pl/Q(T roi \ 
where p b is the 6-type electric dipole moment (much larger than 
the a-type dipole for all molecules considered here). Q(T rol ) is 
the partition function at rotational temperature T rot and is pro- 
portional to (ABC) -0 ' 5 , where A, B, and C are the rotational 
constants of the species concerned. We can compare the inten- 
sities of transitions from the synthesized esters and ethers to 
those in ethanol by assuming that all molecules have the same 
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column and T rol . Then for similar transitions, we have 

TAX) = MX) [ ABC(X) 1° 5 
T^EtOH) /z*(EtOH) LABC(EtOH)J ( } 

Using rotational constants and dipole moments from Hayashi 
& Kawad (1974, 1975, 1980), we find T A (X)/T A ( EtOH) = 0.3, 
0.1, and 0.1 for X — methyl ethyl ether, diethyl ether, and 
methyl propyl ether, respectively. Since our numerical results 
indicate that, for times less than a few times 10 4 yr, the abun- 
dances of large ethers and esters are much less than ethanol, it 
is clear that their detection will prove a formidable challenge to 
radioastronomers. B. E. Turner (private communication) has 
searched unsuccessfully for methyl ethyl ether, dimethyl ether, 
and ethyl formate in Orion-KL, finding upper limits to the 
column densities of about 10 14 cm 2 . In W51 M, he finds the 
upper limit to the diethyl ether abundance is about (2-10)% 
that of ethanol, consistent with the results shown in Figures 3 
and 4 if hot core ages are less than a few times 10 4 yr. More 
sensitive searches for diethyl ether might prove rewarding. 

Interstellar molecules can also be detected through their fun- 
damental vibrations in the infrared, which can appear in 
absorption against a bright background source (i.e., a lumi- 
nous protostar). Saturated aliphatic ethers, for example, have 
two characteristic strong absorption bands in the range 1 1 50— 
1060 cm -1 (8.70-9.45 //m), due to the C— O— C stretching 
vibration (Shimanouchi 1978). The asymmetric mode around 
1 130 cm” 1 (8.85 /mi) is generally the strongest. For pure ethers 
(i.e., diethyl ether), the C”0— C stretching mode occurs at 
about 1 1 10 cm” 1 (9.0 /on) and the asymmetric mode is absent. 
These modes are very strong with an expected integrated inten- 
sity of 3 x 10” 17 cm molecule” 1 . Assuming a thermal line 
width at 300 K, this strength translates into a column density 
of 10 1 2 cm ” 2 in a rotational state for an optical depth of 0.05 in 
one specific rovibrational transition. Hence, for a total absorb- 
ing hot core column density of 10 23 cm” 3 this corresponds to 
an abundance of 10” 11 in one rovibrational state. Emission in 
the rovibrational lines requires high densities and high tem- 
peratures of the gas. The critical density for rovibrational tran- 
sition is on the order of ~10 13 cm” 3 while the excitation 
energy is ~300 K, so that for typical hot core densities and 
temperatures not much collisional excitation is expected. Mol- 
ecules in a dense molecular clump near a luminous protostar 
will resonantly “scatter” the incident IR radiation from the 
dust photosphere. Hence, the molecular rovibrational tran- 
sitions will appear as “emission” lines on top of the IR dust 
continuum from that clump. Typically, we can expect to see an 
optical depth of unity in the dust into the clump. The best 
wavelength region to search is then probably the dust opacity 
gap in the 5-8 fim region. At 8 /zm, outside of the silicate 
feature, this corresponds to a hydrogen column density of 
1.5 x 10 23 cm -2 . Adopting the absorption strength for diethyl 
ether as a typical strength, we find a molecular scattering 
optical depth of 0.5X_ 9 per unit of dust extinction optical 
depth, where X _ 9 is the molecular abundance in the particular 
level under consideration in units of 10“ 9 . The scattered flux 
from a clump at a distance, d , from a protostar with a lumi- 
nosity, L # , and color temperature, T c , is given by (Pendleton, 
Tielens, & Werner 1990) 

(8) 


with ft the solid angle of the clump and x 8ca the scattering 
optical depth. Adopting L*= 10 5 L Q , d = 0.02 pc, (10" at 500 
pc for Orion), T c = 500 K, ft = 3 x 10” 10 sr (2" radius), and 
r sca = 0.5 yields F s = 6 x 10” 16 W cm 2 /mi” 1 . For compari- 
son, assuming the clump is optically thick in the dust and a 
typical temperature of 150 K, the dust continuum flux at 8 /zm 
is calculated tobe7 x 10” 16 W cm” 2 /im _l . Hence, the line- 
to-continuum ratio is about unity. The dust continuum is, of 
course, very sensitive to the adopted dust temperature and 
hence the distance between clump and star. It may actually be 
advantageous to move relatively far from the illuminating star 
where dust temperatures are low and thus much lower molecu- 
lar scattering optical depths could be detected. Absorption or 
emission IR spectroscopy — depending on geometry — may be 
the best way to obtain secure identifications of large inter- 
stellar molecules in hot cores (via their rovibrational spectra). 

5. CONCLUSIONS 

We have demonstrated that evaporation of molecular ice 
mantles containing simple alcohols drives molecular complex- 
ity leading to the formation of pure and mixed ethers, as well as 
some esters. We find that the abundances and relative propor- 
tions of pure and mixed ethers depend upon the initial mantle 
abundance of the precursor alcohol. Ethers formed from 
branched alcohols are more abundant than those formed from 
the unbranched form, assuming that both isomeric forms were 
equally abundant in the ices. We find that CH 3 OC 2 H 5 , 
(C 2 H 5 ) 2 0, and HCOOC 2 H 5 could be abundant in methanol/ 
ethanol-rich hot cores like the Orion Compact Ridge and W51 
M. The abundances of species derived from higher alcohols are 
less certain due to lack of observational constraints, but, in hot 
cores where propanol and butanol are present, i-C 3 H 7 OCH 3 
and (i-C 3 H 7 ) 2 0 will be the most abundant species, followed by 
i-C 3 H 7 OOCH, (n-C 3 H 7 ) 2 0, and r-C 4 H 9 OCH, respectively. 
An important issue is that, as hot core chemistry is of a highly 
transient nature, the observation of an alcohol (e.g., C 2 H 5 OH) 
at a moderate abundance does not exclude it from having been 
much more abundant in the original ices. The presence of 
detectable abundances of CH 3 OC 2 H 5 and (C 2 H 5 ) 2 0 would 
confirm that C 2 H 5 OH can attain large solid-state abundances. 
If the most pessimistic model considered here is the most 
appropriate, then CH 3 OC 2 H 5 will be the only derived mol- 
ecule with an abundance greater than 10“ 1 1 . 

In conclusion, searches for several of the large organic mol- 
ecules considered in this paper should be focused on those 
sources known to be rich in ethanol. Observations are required 
to measure the relative abundances and distributions of meth- 
anol and ethanol in a larger sample of hot core sources. Many 
relevant reactions have been studied in the laboratory and 
several assumptions, concerning rates and products, have 
recently been experimentally confirmed (Smith, Spanel, & 
Millar 1994; Lee et al. 1992). Further laboratory studies are 
needed to determine reaction products and to explore the pos- 
sible importance of alcohol-driven ion-neutral chemistries. 
Theoretical studies of grain surface chemistry, and comparison 
of the results with laboratory ice analog experiments, are 
needed to discover the formation routes to methanol and 
ethanol, as well as to assess if general solid-state pathways to 
interstellar alcohols exist. 
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